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1 Introduction 



Applications of Random Matrix Theory (RMT) to the physics of strong inter- 
actions have a long history (see Chapter 2). RMT was introduced to nuclear 
physics by Wigner to describe the level spacing distribution of nuclei |Wig55| . 
This paper inspired a large body of early work on RMT which is summarized in 
the book by Porter |Por65j . An important conceptual discovery that emerged 
from this field is the large approximation. It first appeared in the work of 
Wigner, and became an integral part of QCD through the seminal work of 't 
Hooft [tHo74j ■ which showed that the limit of a large number of colors is dom- 
inated by planar diagrams. It was realized soon thereafter that combinatorial 
factors can be obtained from matrix integrals |Bre78j . This culminated in the 
random matrix formulation of quantum gravity in 2d, which is a sum over ran- 
dom surfaces that can be triangulated by planar diagrams, (see Chapter 30 
and |DiF93] for a review). 

The main focus of this chapter to apply RMT to spectra of the Dirac op- 
erator both at zero chemical potential, when the Dirac operator is Hermitian, 
and at nonzero chemical potential, when the Dirac operator is nonhermitian. 
Because the Euclidean Dirac operator can be interpreted as a Hamiltonian, this 
application is closer in spirit to the original ideas of Wigner than to the work of 
't Hooft. However, we have benefited greatly from the mathematical techniques 
that were developed for large N QCD and 2d quantum gravity. 

Applications of RMT to QCD have been reviewed extensively in the litera- 
ture |Guh97bl [Ja598bl IVerOOl IVerOSl IAke07j . These papers offer both additional 
details and different points of view. Phenomenological applications of RMT to 
QCD are not discussed in this chapter (see [VerOOl [Ake07| for reviews). 

In the first half of this chapter we introduce chiral random matrix theory and 
its applications to QCD at zero chemical potential. In the second half, chiral 
random matrix theories for QCD at nonzero chemical potential are discussed. 

1.1 Spontaneous Symmetry Breaking in RMT 

One of the essential features of Random Matrix Theory is spontaneous symme- 
try breaking. The real part of the resolvent 

where the average is over the probability distribution of the ensemble N x N 
anti-Hermitian random matrices, D, can be expressed as the replica limit 

ReG(z) = lim -i--^log(Z„(z)), with z £ TZ, (2) 
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where generating function is defined by 



Zn{z) = (det"(L' + z)det"(z - D)). (3) 

For z = 0, the generating function is invariant under Gl{2n). This symmetry 
is broken spontaneously to Gl{n) x Gl{n) by a nonzero value of ReG(z) 

lim lim ReG{z) / 0. (4) 

The order of the limits is essential - the reverse order gives zero. This can be 
seen by expressing the resolvent in terms of eigenvalues of D. 

When we integrate the resolvent over the contour C in the complex plane 
that is the boundary of [— | < Rez < f ] x [— ^Ax < Imz < ^Ax] we obtain 

N j, G{z)dz = 2TTiNAx = 2Trip{0)Ax, (5) 

where iV^x is the number of eigenvalues enclosed by the contour and p(0) is the 
average spectral density around zero. In the limit e — > the l.h.s. is given by 
2iAa;ReG(|). The discontinuity of the resolvent and p(0) are thus related by 

limRe(G(-)) = ^. (6) 

This formula is known as the Banks-Casher formula |Ban80j . It relates the 
order parameter for spontaneous symmetry breaking to the spectrum of the 
associated operator. 

When the spectrum of D is reflection symmetric, i.e. D and —D have the 
same spectrum, then det{z — D) = det(Z? + z), and the generating function for 
the real part of the resolvent is given by 

= (det"(Z5 + z)). (7) 

This is the case when the random matrix ensemble has an involutive symmetry, 
ADA = —D with = 1, which is the case for the Dirac operator in QCD. 

1.2 The QCD Partition Function 

The Euclidean QCD partition function for Nf quarks with mass ruf is given by 

l[{det{D + mf)). (8) 

The average is over gauge fields weighted by the Yang-Mills action. The gauge 
fields are elements of the Lie Algebra SU (Nc) and can be in the fundamental 
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or adjoint representation. The theory that describes the strong interactions 
has Nc = 3 with gauge fields in the fundamental representation. The Dirac 
operator D is a function of the "random" gauge fields. 

Applications of RMT to QCD differ in several respects from other applica- 
tions. First, the physical system itself is already a stochastic ensemble. Second, 
the QCD partition function is not a quenched average, but the fermion determi- 
nant describes the quark degrees of freedom. Third, the average is over Dirac 
operators with different rank. The reason is that, according to the Atiyah- 
Singer index theorem, the number of topological zero modes is equal to the 
topological charge of the gauge field configuration. Because Dirac operators 
with different number of topological zero modes turn out to have spectra with 
different statistical properties, we will treat each topological charge sector sep- 
arately. Fourth, QCD is a quantum field theory that has to be regularized and 
renormalized. We notice that the low-lying Dirac spectrum is gauge invariant 
and renormalizable |Giu09j . 

2 QCD and Chiral Random Matrix Theory 

In this section we construct a RMT with the global symmetries of QCD and 
determine the parameter range for which it is equivalent to the QCD. 

2.1 Symmetries of QCD 

To analyze the global symmetries of QCD we consider the Dirac operator for a 
finite chiral basis. Then it is given by a matrix with the block structure 



where A is a complex x N- matrix. Therefore, all nonzero eigenvalues 
occur in pairs ibAfc. The number of zero eigenvalues is equal to |A^+ — A^_| and 
is interpreted as the topological charge. Paired zeros may occur, but this a set 
of measure zero and of no interest. The block structure of ([9]) is due to the 
axial U{1) symmetry. The partition function also has a vector U{1) symmetry 
related to the conservation of baryon charge. For Nf > 1 the Dirac operator 
is the direct sum of Nf one-flavor Dirac operators so that QCD has the axial 
flavor symmetry UA{Nf) and the vector flavor symmetry Uy{Nf). 

For QCD with Nc'> 2> and gauge fields in the fundamental representation 
there are no other global symmetries. Because SU{2) is pseudoreal, for Nc = 2 
the Dirac operator has an anti-unitary symmetry |Ver94a] 




(9) 



[UK, D]=0 



(10) 
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with K the complex conjugation operator and U a fixed unitary matrix with 
= 1. Then it is always possible to find a basis for which A becomes real 
|Dys62| , and for m = 0, we have that detZ) = det'^^. For A'^^ massless fla- 
vors the quark determinant occurs in the partition function as det^^-'' A. This 
enlarges the flavor symmetry group to U{2Nf). The third case is when gauge 
fields are in the adjoint representation. Then the Dirac operator also has an 
anti- unitary symmetry but now C/^ = — 1 |Ver94aj . In this case it is possible to 
construct a basis for which the Dirac operator can be rearranged in self-dual 
quaternions |Hal95b] . The flavor symmetry is also enlarged to U{2Nf). 

QCD in three Euclidean space-time dimensions does not have an involutive 
symmetry. In that case the flavor symmetry is U{Nf) for gauge fields in the 
fundamental representation. For two colors in the fundamental representation 
and for any number of colors > 2 in the adjoint representation the symmetry 
group is enlarged to 0{2Nf) |Mag99a| and Sp{2Nf) |Mag99b| , respectively. 



2.2 Spontaneous Symmetry Breaking 

The features that mostly determine the physics of QCD at low energy are spon- 
taneous symmetry breaking and confinement. Because of confinement quarks 
and gluons do not appear in the physical spectrum so that QCD at low en- 
ergy is a theory of the weakly interacting Goldstone modes associated with the 
spontaneous breaking of chiral symmetry. 

According to the Vafa-Witten theorem [Vaf83j , global vector symmetries of 
vector-like gauge theories cannot be broken spontaneously. The order parame- 
ter for the breaking of the axial symmetry is the chiral condensate 



lim lim — ( Tr-— 

ma-*0 V— >oo V \ D + iria 



with p{{))/V the spectral density of the Dirac operator per unit volume of space 
time. Because of the Banks-Casher relation, the absolute value of the chiral 
condensate is flavor independent, but its sign is determined by the sign of m^. 
In Table [J we give the symmetry breaking patterns |Pes801 Vys85| for the 



theories mentioned above. We also give the breaking pattern for QCD in 
three dimensions, but refer to the literature for additional discussions |Ver94bt 
Mag99at [Mig99b| EjaOOl |Diiii02| . 



2.3 Chiral Random Matrix Theory 

Since the global symmetries of the Dirac operator are a direct consequence of 
its block structure and the reality properties of its matrix elements, it should 
be clear how to construct a RMT with the same global symmetries: just replace 
the nonzero matrix elements by an ensemble of random numbers. Such chiral 
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Random Matrix Theory (chRMT) is defined by |Shu921 IVer94aj 



z!^M^{{mk}) = I DAlldet ( '^f ) P{A), (12) 

where A is an x matrix, and the integration is over the real (and 
imaginary) parts of A. The reality classes are denoted by the Dyson index (3 
which is equal to the number of degrees of freedom per matrix element and is 
the same as in the corresponding QCD like theory (see Table [1]). The properties 
of the chRMT partition function do not depend on the details of the probability 
distribution. This is known as universality |Bre951 [Jac96al lAke961 IGuh97aj and 
justifies to simply average over a Gaussian distribution 

P(^) = ce-^^'T'-^^'^) with N = N+ + N^. (13) 

Both the Vafa-Witten theorem [Vaf83j and the Banks-Casher formula ([6]) 
apply to chRMT. The global symmetry breaking pattern and the Goldstone 
manifold are therefore the same as in QCD. 

In chRMT, N is interpreted as the volume of space-time. This corresponds 
to units where N/V = 1 so that S can be written as 'SV/N and Eq. (jlSp 
becomes dimensionally correct. Notice that the matrix elements of the Dirac 
operator and its eigenvalues have the dimension of mass. The normalization of 
(|13p is such that S can be interpreted as the chiral condensate that satisfies the 
Banks-Casher relation. 



2.4 Chiral Lagrangian 

The low-energy limit of QCD is given by a chiral Lagrangian and necessarily 
has the name transformation properties as the QCD partition function. The 
Lorentz invariant chiral Lagrangian to 0{M) and 0{p^) is given by 

L = ^F^ Tr df,Udf,U^ - Tr [MU^ + M^U], (14) 

where f/ G G/H with G the global symmetry group that is spontaneously 
broken to H (see Tabled] and F is the pion decay constant. In the domain 

M«^«^ (15) 

the kinetic term factorizes from the chiral Lagrangian [Gas87j . In this domain, 
known as the microscopic domain |Ver94aj . the mass dependence of the QCD 
partition function in the sector of topological charge i' is given by 

Z^{M)= f dC/det'^t/e"^^T'-[^^^'+*^'^]. (16) 
JueG/H 
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Microscopic Domain Cliiral Domain Macroscopic Domain 

chRMT XPT 

Figure 1: Domains for the mass dependence of the QCD partition function. 

The first inequahty in (jlSp can be rewritten as l/M.,^ ^> L, i.e. the pion 
Compton wavelength is much larger than the size of the box. The domain 
where QCD is described by the chiral Lagrangian (jl4p including the kinetic 
term, i.e. where M^r ~ p <C Aqcd; will be called the chiral domain. The zero 
momentum partition function (jl6p and the chiral Lagrangian (jl4p are the first 
term of the e and p expansion |Gas87j . respectively. Therefore these domains 
are also known as the e domain or p domain corresponding to a counting scheme 
where M ~ 1/1/ and M ~ 1/VV, in this order. 

ChRMT can be reformulated identically in terms of a nonlinear cr-model 
even for finite size matrices |Shu921 IJac96aj (see also chapter 7) . In this for- 
mulation the random matrix partition function is given by a integral over two 
types of modes, would be Goldstone modes with mass ~ \/ Nm and massive 
modes with mass ^ ^/N (the number of integration variables does not depend 
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Adjoint 
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0{Nf) X 0{Nf) 
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Table 1: Classification of QCD like theories in three and four dimensions. The 
Dyson index, f3, is the number of degrees of freedom per matrix element {(3 = 1 
and /3 = 4 are interchanged for staggered fermions). 



on N). To leading order in and m the Goldstone mode part factorizes from 
the partition function. Since the pattern of chiral symmetry breaking for QCD 
and chRMT is the same, their mass dependence is the same in the microscopic 
domain (jlSp . The correction terms are of 0{Nm?), so that chRMT results are 
universal for m <^ 1/^/N. With the identification of A'^ as V, this corresponds 
to the e domain. 

In lattice QCD, and QCD in general, the topological charge does not affect 
the block structure of the Dirac matrix but leads to nontrivial correlations 
between the matrix elements. Then it is natural to work at fixed 0-angle (in 
fact at 9 = 0) with partition function given by Z{m,6). However, properties of 
Dirac eigenvalues depend on the topological charge, and in comparing lattice 
QCD and chRMT, Dirac spectra are sorted accordingly. In chRMT, topology 
is included by means of the block structure of the Dirac matrix right from the 
start, and it is natural to work with fixed topological charge. The two partitions 
functions are related by 

oo 

Z{m,0)= J2 e'^^Z-'im). (17) 

u=—oo 

Since the spectrum of the Dirac operator depends on the topological charge it 
makes sense to introduce the topological domain |Leh09j . as the domain where 
the average properties of the eigenvalues are sensitive to the topological charge. 
This domain is expected to coincide with the microscopic domain. 

2.5 Generating Function for the Dirac Spectrum 

The generating function for the Dirac spectrum is also given by d?]) with the 
determinant of the physical quarks contained in the average. Therefore, z 
plays the role of a quark mass and the theory will have Goldstone bosons with 
squared mass 2zTj/F'^. Therefore, for physical quark masses, i.e. quark masses 
that remain fixed in the thermodynamic limit, we can choose |Ver95t IOsb98b) 

^ « ^ ^ ^Th. (18) 

In this domain the z-dependence of the generating function at fixed topological 
charge v is given by the zero momentum partition function (jl6|) with quark 
masses equal to z. The energy scale in (jlSp is known as the Thouless energy. 
A similar conclusion was reached in | Jan98aj . The volume dependence of the 
Thouless energy has been confirmed by lattice simulations |Ber98bl IBer99j . 

The number of eigenvalues that is described by chRMT scales as E'xh/'^A = 
-F^L^/tt. Since F ~ \fNc, this number increases linearly with |Nar04j . 

The determinants containing z have to be quenched which can be done by 
the replica trick (see chapter 8) or the supersymmetric method (see chapter 
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7). A supersymmetric version of the chiral Lagrangian is known and the zero 
momentum integral has been evaluated analytically |Qsb98al IDam98b] . 



2.6 Chiral Random Matrix Theory and the Dirac Spectrum 

In an influential paper that motivated the introduction of chRMT, Leutwyler 
and Smilga |Leu92j proposed to expand the the QCD partition function at fixed 
topology in powers of m and equate the coefficients to the expansion of the same 
ratio for the low energy limit of QCD, i.e. for the zero momentum partition 
function. The sum rules are saturated by eigenvalues in the microscopic domain 
which appear in the combination XkV- 

With eigenvalues that scale as 1/V, the microscopic scaling limit of the 
spectral density can be defined as |Shu92l IVer94a| 

with the microscopic scaling variable defined hy z = WT,. For z ^ ^/V A"^ 
the microscopic spectral density is given by chRMT. For all three values of the 
Dyson index, it can be expressed in terms of the Bessel kernel |For93| 

, X ^ xJa^x{x)Ja{y) - yJa{x)Ja+i{y) 

Ka{x,y) = ^fxy 5 5 . (20) 

The microscopic spectral density for /5 = 2 is given by |Ver93] 

/5f=2'«(x) = hm Kaix.y) = \x{JI(x) - Ja+i{x)Ja-i{x)), (21) 

with a = Nf + The microscopic spectral density for (3=1 and /? = 4 can 
be obtained by rewriting the partition function in terms of skew-orthogonal 
polynomials. For /3 = 1 we find |Ver93l IFor98l lAkeOSbj 



^ ' 4 ^ ' 4 ./n |z — w| \wdw z dz 

= P\ 

with a = 2Nf + and for /3 = 4 the result is |Nag95[ lAkeOSb] 

pf=4-«(^) = 2z^ duu^ f dv{l-v^)v-^/'^Ka{2uz,2uvz) 

f3=2,2a 



Ps 



1 /-akl 

(2z) - -J2a{2z) / dtJ2a{t). (23) 
^ JO 



where a = A^/ + 2|i^|. Similar relations between the microscopic spectral density 
and the kernel for (5 = 2 exist for an arbitrary invariant probability potential 
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Figure 2: The valence quark mass dependence of the chiral condensate. Lat- 
tice data obtained by the Columbia group [Cha95j are compared to analytical 
chRMT result given in Eq. (|25p . 



|Sen981 IKleOOj . and can be exploited to show universality for (3=1 and /? = 4 
from the universality of the microscopic [3 = 2 kernel |Ake96j (see Chapter 6). 

It is of interest to study the critical exponent of the spectral density at a 
critical deformation of the Dirac operator for which /?(0) = which is a different 
universality class. It is unlikely that this critical exponent is equal to the mean 
field value of 1/3 [Jac95j . Other critical exponents can be obtained by fine 
tuning the probability distribution |Ake97l IAke02al IJan02j . 

There are a large number of lattice results for the microscopic Dirac spec- 
trum. The microscopic spectral density (j2ip was first observed for lattice QCD 
Dirac spectra through the mass dependence of the resolvent |Ver95j defind as 
(earlier direct comparisons for the microscopic spectral density were obtained 
for gauge field configurations given by a liquid of instantons and anti-instantons 
|Ver94cp 




The microscopic limit of S(m„), obtained by replacing the sum by an integral 
over the microscopic spectral density, is given by a simple expression in terms 
of Bessel functions |Ver95] 

lNf+\u\{mv)KNf+\y\{m^) + lNf+\u\+i{'>K)KNf+\u\-i{'>K) ■ 

(25) 
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z 



Figure 3: Microscopic spectral density for SU{2) gauge group in the adjoint 
representation (left, taken from |Edw99j ) and in the fundamental representation 
(right, taken from |Ber98aj ). and for SU{'i) gauge group in the fundamental 
representation (middle, taken from |Dam98a] ) . 

In Fig. [2] we compare this result to the chiral condensate (j24p obtained from 
lattice simulations |Cha941 ICha95j for two flavors and various values of the 
coupling constant. The average is over the Yang- Mills action and the fermion 
determinant with the masses in the fermion determinant kept fixed. Agreement 
is found with the quenched result because the physical quark masses are much 
larger than the microscopic scale. There is no dependence on the topological 
charge because the staggered lattice fermions are not close enough to the contin- 
uum limit (This point was investigated in more detail in |Dam99c"| IFar99j . More 
recently it was shown |Won04^ IFol041 IFol05] that dependence on topology as 
predicted by chRMT is reproduced by staggered fermions if we are sufficiently 
close to the continuum limit). Similar lattice results have been obtained for 
(3 = 1 and /3 = 4 |Dam99b] together with the corresponding analytical results. 

The first direct observation of the micropscopic spectral density in QCD was 
made for staggered fermions with two colors in the fundamental representation 
|Ber97aj for which the Dyson index is /3 = 4 (see right panel of Fig. [3|) . Results 
for QCD with A^c = 3 were obtained in |Dam98al IGoc98j (Fig. El middle). 
The left panel of Fig. [3] is for staggered fermions in the adjoint representation 
|Edw99j . 

The distribution of the smallest Dirac eigenvalue is given by P^i^{s) = 
—E'[s) with E{s) defined as the probability that there are no eigenvalues in 
the interval [0, s). In Table [2] we summarize analytical results for the quenched 
case |Ede88[ IFor931 IWil971 IDamOOj . The result for = is particularly simple. 
Expressions for the /c'th smallest eigenvalue at arbitrary quark mass, topological 
charge and Dyson index are known as well |Wil97[ INis98[ IDamOnj . A useful 
measure to compare lattice QCD results with chRMT predictions is the ratio of 
low lying eigenvalues [Giu03t IFod09j . Agreement with the topology and mass 
dependence has become an important tool in lattice QCD to test the lattice 
implementation of chiral symmetry and topology. 
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The joint probability distribution of chRMT only depends on Nf and 
through the combination 2Nf + /3z^. This property, known as flavor-topology 
duality |Ver97j . has been observed in lattice QCD Dirac spectra |Fuk07j . 

Contrary to correlations of low-lying eigenvalues, bulk spectral correlations 
can be investigated by spectral averaging, and for large lattice volumes, the 
Dirac spectrum of a single gauge field configuration is sufficient to obtain sta- 
tistically significant correlators. Excellent agreement with the Wigner-Dyson 
ensembles was obtained |Hal95a] without sign of a Thouless energy. It turns 
out |Guh98] that the Thouless energy scale is due to ensemble averaging. The 
conclusions is that there is no spectral ergodicity beyond the Thouless energy. 





1/ = 


V = 1 


general v 

V odd for (5=1 


p = l 

[5 = 2 
[5 = A 


C 

i(eC(^_l)+e-C(C+l)) 




|det/i_j+2(C) 



Table 2: Results for the distribution of the smallest Dirac eigenvalue for the 
quenched case. In the last row Oj is an expression in terms of a sum over 
partitions of j. For explicit expressions we refer to [Ber98a| . 

2.7 Integrability 

For [5 = 2 the partition function of invariant random matrix theories can be 
interpreted as a partition function of noninteracting fermions which is an inte- 
grable system. This is the reason that the low energy QCD partition function, 
the unitary matrix integral (jl6p . obeys a large number of remarkable relations. 
The Nf flavor partition function in the sector of topological charge v can be 
written as [BroSlal IBroSlbl IGuh96l IJac96bj 



In the limit xj. 



det \x' 



k-iAi-i) 
k 



Xk = nikVT.. (26) 



A({xi}) 

X this partition function reduces to a Hankel determinant 

k+l ■ 



Z]^^, =det[(x9,)'^+7,(x)]fc,^=o,...,^,-i. (27) 

Applying the Sylvester identity |For02| relating the determinant of a matrix to 
co-factors gives the Toda lattice equation |Kan021 Spl03a| 



{xd^f\ogZ'{,(x) = 2Nfx' 



z- 



(28) 



After taking the replica limit of this recursion relation we arrive at the following 
compact expression for the resolvent |Spl03a| 



xdxxG{x) = lim log Z^r (x) 



2x^Z^{x)Zl^{x). 



(29) 
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This factorized form is a general property of the spectral density and correlation 
functions of RMT's with (3 = 2. Application of the replica limit to a discrete 
recursion relation does not require analyticity in the replica variable and this 
way problems with the replica limit can be circumvented |Ver851 IKan02j (see 
Chapter 8). 

As a consequence of integrability relations, the zero momentum partition 
function satisfies Virasoro constraints. They provide efficient way to determine 
the coefficients of the small mass expansion |Dam99at IDalOl] . 

In addition to the relations discussed above we would like to mention the 
following relations: i) The Toda lattice equation can be formulated for finite 
size random matrices by exploiting the properties of orthogonal polynomials 
|Ake04] . ii) In the microscopic domain the partition functions in 3 and 4 di- 
mensions are related by |Ake99[ lAkeOnbi IAndn4j 

ZQc/D,iM) = Z''=-'/\{x,})Z''='/\{x,}). (30) 

iii) A;-point spectral correlation functions can be expressed into partition func- 
tions with pk additional flavors |Ake98j . iv) The correlation functions of in- 
variant RMTs can be expressed in terms of the two-point kernel. This leads 
to consistency relations between various partition functions |Ake98j v) Inter- 
preting the quark mass as an additional eigenvalue leads to relations between 
correlators of massive and massless partition functions [AkeOOa] . 

3 ChRMT at Nonzero Chemical Potential 

An important application of chRMT is to QCD at nonzero chemical potential 
fj,. In that case the Dirac operator is given by 

Dill) = Dili = 0) + ^70. (31) 

Since D{ii = 0) is anithermitian, D[ii ^ 0) has no hermiticity properties, and 
its eigenvalues are scattered in the complex plane. Because the determinant 
of the Dirac operator is complex, the QCD partition function at ^ 7^ is the 
average of a complex weight, and unless the chemical potential is small, it can- 
not be simulated by Monte-Carlo methods. For that reason chRMT has been 
particular helpful to answer questions that could not be addressed otherwise. 
In particular, the following issues have been clarified: i) The nature of the 
quenched approximation |Ste96b] . ii) The relation between the chiral conden- 
sate and the spectrum of the Dirac operator for QCD with dynamical quarks 
[OsbOSj . iii) The expectation value of the phase of the fermion determinant 
|Spl06| . iv) The low-energy limit of phase-quenched QCD and the spectrum 
of the Dirac operator [TouOOj . v) The geometry of the support of the Dirac 
spectrum [TouOOl lOsbOSb] . 
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The QCD partition function at nonzero chemical potential is given by 

ZQCD = {Yldet{D + mk + fiklo))- (32) 

k=l 

where the average is over the Yang-Mills action. The chemical potential for 
different flavors is general different. Two important special cases are: /Xfc = /i 
when// is the baryon chemical potential, and the case for an even number of 
flavors Nf = 2n with fi/^ = fi for k = 1, - ■ ■ ,n and = —fj, for k = n+1, • • • , 2n. 
In the second case the partition function is positive definite because 

det{D + m- fi-fo) = det{D^ + m + fijo) = det*{D + m + /i7o). (33) 

For n = 1, /_f can be interpreted as an isospin chemical potential [SonOOj . Since 
the determinant appears together with its complex conjugate, this partition 
function is also known as the phase quenched two-flavor partition function. 

3.1 Dirac Spectrum 

A particular useful tool for studying the spectrum of a nonhermitian operator 
is the resolvent ([1]). The spectral density is given by 

p{z, z*) = 6\z - A,)) = (34) 

k 

and can be interpreted as the two-dimensional electric field at z of charges 
located at A^. When the eigenvalues are on the imaginary axis this picture 
illustrates that G{z) has a discontinuity when z crosses the imaginary axis. 
When eigenvalues are not constrained by Hermiticity, because of level repulsion, 
they will scatter into the complex plane. Using the electrostatic analogy, the 
resolvent will be continuous. If z is outside the spectrum, G{z) is analytic in z. 
The resolvent cannot be expressed in terms of Eq. ([7]) but rather as 

G{z) = Yim-^Zn{z,z*), (35) 
n^o ndz 

where Z„(z,z*) is the phase quenched partition function |Ste96bj 

Zn{z, z*) = (det"(L' + z)dei'\D^ + z*)). (36) 

Lattice QCD Dirac spectra were first calculated in |Bar86j . As remarkable 
features we note that the spectrum is approximately homogeneous, and that it 
has a sharp edge which both are explained by chRMT. 
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3.2 Low-Energy Limit of QCD and Phase Quenched QCD 



According to the definition of the grand canonical partition function, the free 
energy at low temperature does not depend on the chemical potential until it is 
equal to the lightest physical excitation (per unit charge) with charge conjugate 
to fj,. For QCD this implies that the chiral condensate at zero temperature 
does not depend on // until n = mj\f/3 (with uin the nucleon mass). The Dirac 
spectrum, is n dependent, though, which seems to violate the Banks-Casher 
relation |Ban80j . This problem is known as the 'Silver Blaze Problem' [Coh03j . 

At nonzero isospin chemical, fij, the critical chemical potential is equal to 
fij = mTt/2. Beyond this point, pions will Bose condense. For light quarks, this 
phase transition can be studied by chiral perturbation theory, and for quark 
masses in the microscopic domain it is described by chRMT. At nonzero fij, 
the 'Silver Blaze Problem' is that at zero temperature the chiral condensate 
remains constant until fii = while the spectral density depends on ///. 

The solution is easy: according to the electrostatic analogy, the 'electric field', 
i.e. the chiral condensate, is constant outside a homogeneously charged strip. 
This implies that the width of the strip is determined by the relation fij = 
|Gib861 [TouOOj . Indeed, in terms of eigenvalues, the critical point is when the 
quark mass hits the boundary of the spectrum. 

3.3 Chiral Lagrangian at Nonzero Chemical Potential 

Chiral symmetry remains broken at small nonzero chemical potential. There- 
fore, also in this case, the low-energy limit of QCD is given by a theory of 
weakly interacting Goldstone bosons. As is the case at zero chemical potential, 
the UiiNf) X UR{Nf) invariance of the partition function is broken sponta- 
neously to Uy{Nf). The invariance properties of QCD should also hold for 
the Lagrangian that describes the low-energy limit of QCD. In particular, be- 
cause the chemical potential is an external vector potential, it only enters in 
the combination of the covariant derivative |Kog99| 



Together with the mass term, the 0{p^) chiral Lagrangian is thus given by 



Eq. (|38p shows that the chiral Lagrangian is determined by two constants. 
Since the Dirac spectrum at ;U 7^ is also determined by two constants, the 
eigenvalue density and the width of the spectrum, we can extract the low energy 
constants from the geometry of Dirac spectrum. 

It has been argued that at sufficient large chemical potential QCD will be 
in a color-flavor locked phase with spontaneously broken color-flavor symmetry. 



with B„ = diag({;Ufc})(5i.,o- 



(37) 




(38) 
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In this phase the chiral condensate vanishes, and next order terms in the chiral 
expansion, which are quadratic in the quark mass, have to be taken into account. 
Universal results are obtained by scaling the Dirac eigenvalues with \/V and 
Leutwyler-Smilga sum rules have been derived both for QCD with three colors 
|Yam09j and QCD with two colors |Kan09j . 

3.4 Chiral Random Matrix Theories at /^i 7^ 

In a suitably normalized chiral basis the Dirac operator at nonzero chemical 
potential has the block structure 

id + ji 



^("' = 1 I . (39) 



A chiral random model at nonzero chemical potential is obtained |Ste96bj by 
replacing the matrix elements of d and by an ensemble of random numbers 
exactly as in section 12.31 Also in this case, because of expected universality 
[Ake02b] . it is justified to simplify the model by choosing a Gaussian distribu- 
tion. 

As is the case for /u = 0, we can distinguish three different nonhermitian 
chRMTs jHal97j . with complex matrix elements (/3 = 2), with real matrix 
elements (/3 = 1), and with self-dual quaternion matrix elements (/? = 4). 
They apply to the same cases as discussed in table [TJ The full classification 
of nonhermitian ensembles is based on the Cartan classification of symmetric 
spaces (Zi?96l IBernil |Mag07 . 



The random matrix model (j39p is not unique. Adding /i in a different way 
results in the same chiral Lagrangian as long as the invariance properties of the 
matrix model remain the same. One drawback of the model (I39p is that the 
overall unitary invariance has been lost so that methods that rely on the joint 
probability distribution of eigenvalues cannot be used. A model that does have 
a representation in terms of eigenvalues is defined by |Osb04| 



where C and d are complex random matrices with the same distribution. 
A rerun of the arguments of |Gas871 IDam06j in the microscopic domain 

mV < 1 and <C 1, (41) 

shows that the partition function corresponding to the chiral Lagrangian (I38p 
factorizes into a zero momentum part and a nonzero momentum part. The 
chemical potential and mass dependence reside in the zero momentum part. 

In the microscopic domain, invariance properties of QCD at 7^ that rely 
on global symmetries can also hold for chRMT at /i 7^ and give the same 
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invariant terms in the zero momentum sector. Therefore, in this domain, the 
QCD partition function is given by chRMT at fi ^ 0. Mean field studies only 
involve the zero momentum part of the chiral Lagrangian. Therefore, mean 
field results |KogOO| can also be derived from chRMT. 

Applying the above arguments to the generating function for the Dirac 
spectrum we obtain the zero momentum partition function 



Z';l{z,z*;n) = c [ dUdef'Ue 



UeU{2n) 

with B = S3 and = ( g z* ) ' ^^^^ 

At the mean field level the resolvent is independent of the replica index and 
the partition function can be analyzed for re = 1. For the resolvent we find, 

G(.) = S and Giz) = ?^^±fi (43) 

for Re z > 2/x^F^/S and Re z < 2/x^F^/S, respectively. The eigenvalues are 
therefore distributed homogeneously inside a strip with width 4F^//^/S. In 
agreement with lattice simulations [Bar86j . the eigenvalue density has a sharp 
edge whereas the resolvent is continuous at this point. If z and z* are interpreted 
as quark masses, the squared mass of the corresponding Goldstone bosons is 
equal to = (z + z*)S/F^. The condition Rez < 2^^F^/S can then be 
written as fi = ma/ 2, in agreement with physical considerations. 

3.5 Integrability of the Partition Function 

Remarkably, as was the case = 0, the zero momentum partition function (j42p 
can be rewritten in terms of a Hankel like determinant |Spl03a Spl03b| 



Z;^(z,z*,/i) = n„(zz*)"(i-")det[(za,)'=(z*a,0'^r(^,^*,/^)]M=o,i,-,n-i. (44) 



This form responsible for the integrable structure of the partition function (j42p . 
Most notably, it satisfies the Toda lattice equation 

zB z*3 , loe Z'^(z z* u) - '''' .^*.2 Z^+i(^, z* , ^i)Zi:_,iz, z* , ^i) 
zd^z d..logZ„(z,z ,/.) - —{zz ) [Zii(z,z*,f,)]^ ■ 



which is obtained by applying the Sylvester identity to the determinant in (j44p . 
This equation can be extended to imaginary chemical potential and the two- 
point correlation function |Spl03b , IDamOSl IDamOG] . For imaginary chemical 



potential there is no transition to a Bose condensed state and the re-dependent 
part of the free energy vanishes after differentiation with respect to the masses. 
The nontrivial result in the free energy is O(re^) which gives the two-point 
correlation function. 
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Figure 4: Left: Radial microscopic spectral density for quenched QCD at /i 7^ 
|Wet04j . Right: Spectral density of the overlap Dirac operator as a function of 
the distance, x, to the imaginary axis |Blo06j . 



The replica limit of the Toda lattice equation results in the spectral density 

p{z,z*,n) = -^-^-^ log Z!^{z,z*,fj,) = ^zz* Z'({z, z* , fi)Z'i^{z, z* , fi), (46) 

which was derived in |Spl03b| . The fermionic partition function can be obtained 
by an explicit evaluation of the integral over U{2). The result is given by 

= -e^^^V /"'dAAe-2^^VA^j^(A^Sy)/,(Az*Sy). (47) 

The evaluation of the bosonic partition is more complicated. The inverse com- 
plex conjugated determinants can only be represented as a Gaussian integral 
after combining them into a Hermitian matrix. In order to obtain a convergent 
integral, the Hermitian matrix has to regularized by a mass ~ e. This proce- 
dure is known as Hermitization |Jan971 IFei97j . It turns out that the partition 
function is logarithmically divergent in e. This divergence is due to a single 
eigenvalue close to z and is present even if z is outside the support of the Dirac 
spectrum |Spl03b Spl08| . Because of the Vandermonde determinant, the prob- 
ability of finding two eigenvalues close to z does not diverge. The partition 
function for one pair of conjugate quarks can also be written as an integral 
over Goldstone bosons. Instead of an integral over U{2), using an extension of 
the Ingham-Siegel integral |Fyo01| , we obtain an integral over the noncompact 
manifold of positive definite Hermitian matrices Q 

Z\iz,z*;i,) = hmC, / ^^(Q)eT4^¥C"(Q-^Q-^^)-^^V[Q,<.3][Q-,<x3]] ^ 

w.h..(_»^j), c^(;. :). (48, 

The integral over Q can be performed analytically resulting in 

= 'v^V 4l.%t'^ (49) 



19 



In Fig. m we compare the expression for the spectral density to quenched 
lattice simulations |Wet04j . The lattice data in the left panel of Fig. [Hare in 
the strong-nonhermiticity domain where Re(z)5]/2/i^i^^ < 1 and fi'^F'^V ^ 1. 
In this domain the analytical result can be simplified to |Spl03b , IVer05| 

2 

PNf=o,u{z,z*,u) = —u^zz*Ky{uz*z)Iu{uz*z), (50) 

J TT 



with u defined as 



4^2 ^2 • 



(51) 



The topological index of the staggered lattice Dirac operator at relatively strong 
coupling in the left panel of Fig. [J] is zero. Lattice results for nonzero topolog- 
ical charge have been obtained using the Bloch-Wettig overlap Dirac operator 
|Blo06) and are compared to the analytical expression (j46p in the right panel 
of Fig. H 

Using superbosonization [Hac951 IBun07l IBasOTj , the fermionic and bosonic 
partition function can be combined into a super symmetric partition function 
|Bas07j which can be used to derive the low-energy limit of the generating 
function of the QCD Dirac spectrum at /u 7^ 0. 

3.6 Spectral Density at /i 7^ for QCD with Dynamical Quarks 

Although the spectral density of the Dirac operator for QCD with dynamical 
quarks at // 7^ was first derived using complex orthogonal polynomials |Osb04j . 
a simpler expression is obtained from the Toda lattice equation |Ake04j . 

u t * ^ *T\( 2 2^ Z-=-Hz, z*,{mf},f,)Z^^=^'^f jz, z* , {m^}, /.) 
PnAz,z zz [[{mf-z) '—j^ . 

In Fig. [5] we show a 3d plot of its real part. In addition to a flat region there is 
a strongly oscillating region with oscillations with an amplitude that increase 
exponentially with the volume and period that goes like ~ 1/V. This region 
is absent in the quenched or phase-quenched case and is responsible for the 
discontinuity in the chiral condensate. As we can see from Fig. \5\ we can 
distinguish three phases in the Dirac spectrum. A phase with no eigenvalues, a 
phase with a constant eigenvalue density, and a phase with a strongly oscillating 
eigenvalue density. These phases can be obtained |Osb08bj by means of a mean 
field study of a partition function with masses m, z and z* , similar to the 
analysis of 3 fiavor QCD at nonzero chemical potentials [KogOl . The phase 
diagram is shown in the right panel of in Fig. [5l 

As we argued before, in the microscopic domain, the chiral condensate does 
not depend on the chemical potential. Both the fiat region and oscillating 
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Dirac spectrum for Full QCD. 

Figure 5: Left: The real part of the spectral density of the QCD Dirac operator 
for one-flavor QCD at ^ 0. For better illustration the z-axis has been clipped. 
Right: The phase diagram of the Dirac spectrum 

region give rise to a /i dependent contribution to the chiral condensate, but 
the //-dependence cancels in their sum |Osb05| . This solves the "Silver Blaze 
Problem" [Coh03j and can be explained [OsbOSaj in terms of orthogonality 
relations of the complex orthogonal polynomials. 

3.7 The Phase of the Fermion Determinant 

ChRMT can be used to study the complex phase of the fermion determinant. 
The average phase factor may be calculated with respect to the quenched, 
the phase quenched, or the two-flavor partition function. The phase quenched 
average is given by 



/ ^ / det(D + m + ^i-fp) \ _ ^i+i , , 

^'^''-\det{D^ + m + f,jo)/,^,.~ Z,^,,- ^'^^ 

Since ^i+i does not depend on /x, the latter ratio follows immediately from the 
expression from Z\j^\* given in (j42p . 

The quenched average phase factor can be re-written in terms of a determi- 
nant of complex orthogonal polynomials. In the microscopic domain, the result 
is the sum of a polynomial in ^ and a part with an essential singularity at 
/i = 0. Therefore it cannot be obtained from analytical continuation from an 
imaginary chemical potential which is polynomial in li} [DamOS^ Spl06 lBlo08| . 



3.8 QCD at Imaginary Chemical Potential 

Random matrix models at imaginary chemical potential are obtained by re- 
placing /i — > z/i in the Dirac operator. Then the Dirac operator becomes anti- 
hermitian with all eigenvalues on the imaginary axis. There are two ways of 
introducing an imaginary chemical potential, either as a multiple of the identity 
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or as a multiple of a complex random matrix ensemble. Both models have the 
same symmetry properties and lead to the same universal partition function 
in the microscopic domain. Spectral correlation functions can be obtained by 
means of the Toda lattice equation |Dam06j . or in the second model, by means 
of the method of bi-orthogonal polynomials |Ake08aj . 

Parametric correlations of Dirac spectra in the microscopic domain depend 
on two low-energy constants, F and S, which can be extracted from correlations 
of lattice QCD Dirac spectra |Dam06t |Deg07[ lAkeOSa] . 



4 Applications to Gauge Degrees of Freedom 



The Eguchi-Kawai model Egu82 is the lattice Yang-Mills partition function 
with all links in the same direction identified. In the large limit this model 
is an integral over C/(A^c)-matrices. Although the original hope, that Wilson 
loops of Yang-Mills theory are given by this reduced theory is incorrect, the 
model continues to attract a considerable amount of attention. 

For d = 4 the Eguchi-Kawai model cannot be solved analytically, but for 
(i = 2 it is known as the Brezin-Gross-Witten model |Bre801 IGroSOj 

Z=[ dC/e^^^(^+^^ (54) 

and is identical to the zero momentum partition function (jl6p (see Chapter 
17 for a dsicussion of such group integrals). In the large Nc limit this model 
undergoes a third order phase transition at g'^Nc = 2. 

In the large Nc limit eigenvalues of Wilson loops can be analyzed by means 
of RMT methods. It was shown in [DurSOj that Wilson loops in two dimen- 
sions undergo a phase transition for A'^c — > c« at a critical value of the length 
of the loop. In one phase, the eigenvalues of the Wilson loop are distributed 
homogeneously over the unit circle, whereas in the other phase, they are local- 
ized at zero. This phase transition has been observed in lattice QCD simula- 
tions |Nar06j and has been analyzed in terms of shock solutions of the Burgers 
equation [BlaOSl IBla091 INeuOSj . It can be studied by analyzing the eigenvalue 
distribution of products of unitary matrices |Gud03|. ILohOSj . 



5 Concluding Remarks 

Random Matrix Theory has changed our perspective of the QCD Dirac spec- 
trum. Before the advent of chiral Random Matrix Theory, the discrete structure 
of the Dirac spectrum was viewed as random noise that will go away in the con- 
tinuum limit. Now we know that Dirac eigenvalues show intricate correlations 
that are determined by chiral random matrix theory with one or two low-energy 
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constants as parameters. This implies that we can extract the chiral condensate 
and the pion decay constant from the distribution of individual eigenvalues. 

Chiral Random Matrix Theory primarily applies to the Dirac spectrum, 
and therefore we have to distinguish QCD at zero chemical potential, when the 
Dirac operator is anti-Hermitian, and QCD at nonzero chemical potential with 
a nonhermitian Dirac operator. In the Hermitian case the statistical proper- 
ties of the low-lying Dirac eigenvalues are completely determined by the chiral 
condensate. In the nonhermitian case they are determined by two parameters, 
the chiral condensate, and the pion decay constant. Since the nonhermitian 
Dirac spectrum has the geometry of a strip, the two low-energy constants are 
determined by the eigenvalue density and the width of the spectrum. 

For imaginary chemical potential, the Dirac operator is Hermitian. Al- 
though spectral correlations are completely determined by the chiral conden- 
sate, this is not the case for parametric correlations, which are correlations of 
eigenvalues for two different values of the chemical potential. They require both 
the chiral condensate and the pion decay constant as input parameters. 

Chiral Random Matrix Theory applies to the low-lying Dirac spectrum. 
The scale is set by the momentum dependent terms in the chiral Lagrangian. 
Physically, it is the scale of the quark mass for which the Compton wave length 
of the Goldstone bosons is much larger the size of the box. The scale of the 
chemical potential should also be well below the inverse box size. 

There are two mechanisms to explain confinement in QCD: by condensation 
of monopoles, or by the disorder of gauge fields. The success of Random Matrix 
Theory points to the second mechanism. It our hope that the work discussed 
in this chapter will contribute to the solution of this problem. 
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